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ABSTRACT 

M85OT2006-1 is the latest and most brilliant addition to the small group 
of known Luminous Red Novae (LRNe). An identifying characteristic of the 
previously detected events (M31 RV, V4332 Sgr & V838 Mon) was a spectral 
red-ward evolution connected with an emerging infrared component following the 
optical decay. Here we report on the discovery of a similar feature in Keck/NIRC 
and Spitzer photometry of M85 OT2006-1 six months post-eruption. We find 
that its 2.1—22 /iin spectral energy distribution is best described by a black body 
with effective temperature T e g = 950 ± 150 K and bolometric luminosity L = 
2.9lo;j5 x 10 5 L Q . Assuming spherical geometry, the black body effective radius, 
R = 2.01q;4 x 10 4 R©, and corresponding expansion velocity, v = 870l 2 Jokm s _1 , 
are remarkably similar to the properties of M31 RV 70 days after its eruption. 
Furthermore, we propose a search strategy for LRNe in the local Universe making 
use of the longevity of their infrared excess emission and discuss the expected 
number of events in the Spitzer Infrared Nearby Galaxies Survey. 

Subject headings: stars: individual (M850T 2006-1, M31 RV, V838 Mon), stars: 
variables: other 



1. Introduction 

M85 OT2006-1 was discovered by the Lick Observatory Supernova Search on 2006 Jan- 
uary 7. It appeared as a new transient projected in the outskirts of the lenticular Virgo 
cluster galaxy M85. The outburst, as observed in the optical, exhibited a plateau with 
a duration of about two months and reached an absolute peak magnitude of My — — 13 



-2- 



( IKuIkarni et al 20061 ). Spectroscopic observations indicated that this was neither a classical 
nova nor a sub-luminous supernova. Instead, its early-time spectrum was consistent with a 
black body of T c g ~ 4600 with narrow Ha and Hf3 emission. Following the plateau-phase, 
the optical brightness decreased rapidly and the peak emission evolved towards longer wave- 
lengths. 



A similar evolution has been observed in three ot 



rer transients so far; M31RV in An- 



sively studied V838 Mon ( 



dromeda ({Rich et al.1989) and probably V4332 Sgr (IMartini et al. 19991) and the exten 



Brown et al. 20021 ; iMunari et al. 2002al : IrCimeswenger et al. 2002 



Bond et al. 20031 : lEvans et al. 2003h in the Milky Way. Among this small group, M85 OT2006- 



1 stands out as the most distant and brightest event approaching a bolometric peak lumi- 
nosity of about 5 x 10 6 L Q . 

While the peak luminosities of the four sources differ, their overall observed properties 
are indicative of a common eruption mechanism. Thus, they were suggested to be the first 
members of an emer ging class of eruptiy e transients which was recently dubbed Luminous 
Red Nova^ (LRNe; Kulkarni et al 2006 ). The cause of the peculiar outbursts still needs 
to be fully established, although various explanations were proposed. The most favorable 
scenarios appear to be a stellar merger event in a binary system w ith small secondary-to- 
primary mass ratio (ISoker &: Tylenda 20031 ; ISoker fc Tylenda 20061 ) or a planetary capture 
dRetter fc Marom 2003^ 



One of the most remarkable features of M31 RV and V838 Mon was a strong infrared ex - 
cess which developed a few months after the eruption ( Mould et al. 1990 ; Lynch et al. 2004 1^1. 
This excess emission is commonly associated with condensation of newly formed dust in the 
expelled envelope around the progenitor star. In this paper we present the results of a near 
and mid-infrared search for a similar component in M85OT2006-1. 



2. Observations 



We observed M85O T2006-1 using the Spitzer Space Telescope Infrared Array Camera 
(IRAC; iFazio et al. 20041 ) a t 3.6, 4.5, 5.8 and 8 /im and the peak-up imaging mode of the 
Infrared Spectrograph (IRS; iHouck et al. 20041 ) at 15.8 and 22 /im in July 2006. See Tabled] 



1 The terms Lumin ous Red Variables ( Bryan fc Rover 1992 ). stars erupting into cool supergiants (SECS) 
( Munari et al. 2002rj ) and Mergebursts ( Soker fc Tylenda 20061 ) were suggested as well. 



2 No infrared observations were obtained for V4332 Sgr within the first four years. It is thus unconstrained 
whether it showed a similar evolution or not. 
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for a summary of observations. Standard pipeline post-Basic-Calibrated-Data products were 
obtained from the Spitzer archive. 

A color-composite of the 3.6, 5.8 and 8 /im emission is shown in Figured] (left). Prior to 
undertaking photometry we subtracted the local host galaxy emission estimated using the 
IRAF task mkskycor (Figure [1] right). The flux densities are measured through l'/9 and 3" 6 
apertures for IRAC and IRS, respectively. Aperture corrections were applied according to 
documented values q See Tabled] for a summary of flux densities and apparent magnitudes. 

In addition, we obta ined simultaneous iY-ban d imaging with the Keck Observatory Near 



Infrared Camera (NIRCi lMatthews fe Soifer 1994 ) under excellent (FWHM=0 / /4) conditions 



The data were reduced using the near-infrared processing package IRtools for IRAF (D. 
Thompson, private communication). The host galaxy emission was subtracted as described 
above and the flux c alibration was derived using observations of the standard star SJ 9145 



( jPersson et al. 19981 ). The flux measurement was obtained using a circular aperture of (X'8 



radius. 



3. Results 
3.1. fT-band light curve 



In Figure[2]we show the i\~-band light curve combining the photometry from (IKulkarni et al 20061 ) 
with the NIRC observation presented here. The sparse data are consistent with a linear decay 
at a rate of ~ 0.02 mag per day followed by a flattening at about 120— 180 days. 



3.2. Spectra Energy Distribution 

The 2.1 to 22/im photometry at ~180 days is consistent with a black body with effective 
temperature, T eff = 1030 ± 50 K {x 2 /dof = 3.3/5; solid line in Figure ED- A fit to the IRAC 
and IRS data alone suggests a slightly lower temperature, T eff = 9001*00 K (y 2 /dof = 0.9/4; 
dashed line in Figure [3]) and an additional, hotter component leading to a flux surplus at 
2.1 /im. For the following discussion we will assume a black body temperature of T c g = 
950 ± 150 K for the infrared excess emission. 

The bolometric luminosity after 180 days (as traced by crT^ s ) was L = 2.9^5 x 10 5 L Q 



3 IRAC: 



http://swire.ipac.caltech.edu/swire/astronomers/publications/SWIRE2jdoc_083105.pdf 



http: / / ssc.spitzer.caltech.edu/irs / pu_fluxcal.txt 



IRS: 
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(see Table |2]). The corresponding black body radius was R = [L/(47r<TsT e 4 ff )] 1 / 2 = 2.0^^4 x 
1O 4 R . Using 2006 January 7 as the onset of the eruption, we estimate a mean expansion 
velocity of 870^gQkm s _1 . The total emitted radiation energy is ~ 10 47 erg, dominated by 
the optical plateau during the first three months. 



Discussion and Conclusion 



We have presented the discovery of a strong 3.6-22 /zm excess in M85OT2006-1 at 
~180 days. This thermal infrared component s uggests dust conde nsation in the matter ex- 
pelled during the eruption, similar to M31 RV ([Mould et al. 19901 ) and V838 Mon 
( iKimeswenger et al. 20021 ; iLynch et al. 2004 ) . 



The derived photospheric properties of M85OT2006-1 are resembled closest by those 
of M31 RV. For the latter a 1000 K dust shell with a radius of ~ 8000 R Q was reported 



after 70 days (IMould et al. 1990L see Table [2] for a comparison). Similar to M85OT2006-1, 
M31 RV showed a nearly perfect black body component with only a faint excess at shorter 
wavelengths. Also V838 Mon disp layed a single-component black body emission during the 
first 120 days (cf., iTylenda 20051 ). At this time it's effective temperature, T e s = 3000 K, 
and radius, R eS = 2500 R , differed significantly from those of M85OT2006-1 and M31 RV, 
though. However , V838 Mon was init ially hotter (T e g- ~ 7000 K) and exhibited at least two 
eruption phases (IMunari et al. 2002a ). Furthermore it resides possibly in a binary with a 
bright B3 V star (IMunari et al. 2002a ). making the interpretation of the observed quantities 
difficult. 

Assuming a spherical geometry for M85OT2006-1, we inferred a mean expansion ve- 
locity of ~ 900 km s _1 in June 2006. This value is surprisingly similar to the one de- 
rived for M31 RV at comparable T e g (~ 920km s -1 , see Table |2]). However, 900km s -1 in 
M85OT2006-1 is significantly lar ger than the velocity obtained from the Ha line width in 



February 2006 (350 ± 140 km s ; iKulkarni et al 20061 ). A solution to this apparent discrep- 



ancy lies in the source geometry. In case of an a-spherical expansion, e.g. an oblate geometry 
observed face-on, the velocity (oc black body radius), may be overestimated. Alternatively, 
the large velocity inferred from the infrared observations, can be caused by radiative acceler- 
ation as a result of a long lasting engine activity or multiple phases of energy injection into 
the expanding material. Furthermore, the Balmer emission may have originated in a low 
velocity region unrelated to the matter which produced the late time black body component. 



In V838 Mon the shell became partly transparent after ~120 days. Furthermore, the 
temperature slowly increased again after ~250 days, possibly as the result of a gravitation- 
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ally induced collapse of the inflated envelope (ITylenda 20051). A similar beha vior was also 
reported for V4332 Sgr about nine years after its eruption (ITylenda et al. 20051 ). Whether or 
not M85 OT2006-1 will evolve alike, needs to be addressed with future infrared observations. 

Most of our knowledge on LRNe prior to our analysis was based on detailed studies 
of a single event, V838 Mon. Here we have shown that variations in the infrared evolution 
exists among the small sample of sources. Whether the cause of these variations lies in 
the progenitor, environment or some other unnamed parameter remains to be solved. In- 
deed, already the unde rlying stellar popula tions are not unif orm. While V838 M on resides 
within a B-star c luster (Ufsar fc Bond 2006^. M85OT2 006-1 fcurkarni et al 2006h . M31 RV 
( IRich et al. 19891 ) and V4332 Sgr (ITylenda et al. 20051 ) probably originated from low-mass 
stars. However, it is clearly premature to speculate on a possible bimodality here. 

The next step forward requires a significant increase of the LRNe sample. Estimates 
based o n the few known even ts suggest rates of approx. a dozen per year out to a distance of 
20Mpc (IKulkarni et al 20061 ). Simila r predictions come from the theory of violent mergers 
of both massive and low mass stars (ISoker fc Tylenda 20061 ). Thus, a dedicated search for 
LRNe is expected to net a substantial number of new events and may revolutionize our 
understanding of these elusive transients. 

A promising strategy to detect new LRNe similar to M85 OT2006-1 and M31 RV is 
to make use of their characteristic infrared emission properties. The IRAC colors of a 
900 — 1000 K black body can easily be discriminated from the majority of known stellar 
objects, AGN and unresolved galaxies (Figure 0] left). Only early-type T-dwarfs (T2-T4; 
Patten et al. 20061 ). which have similar T e g, exhibit comparable 3.6— 8 /zm colors. However, 
the latter can be identified by their stron g molecular absor ption bands in the near-infrared 
(e.g., CH 4 , H 2 0; for a recent review see iKirkpatrick 20051 ). which also lead to significant 
deviations from a perfect black body spectrum (offset from black body line in Figure H])Q 

The main advantage of searching in the infrared is that the detection probability is in- 
creased with respect to optical searches by the ratio of infrared-brig ht time (a few years) to 
optical plateau duration (a few months). Using the predictions of (ISoker & Tylenda 20061 ) 
of approx. one LRN per galaxy every 10-50 years, we can estimate the number of antici- 
pated events in an existing sample like the Spitzer Infrared Nearby Galaxies Survey (SINGS; 



4 We note that V838 Mon revealed strong absorption bands (e.g., H 2 0, CO, OH, SH: iLvnch et al. 2004h 
after ~ lyear together with a hot continuum (Toff ~ 2700 K) inconsistent with a simple black body. Thus, 
it would have been difficult to identify by its infrared colors alone. Whether this is caused by intrinsic or 
environmental differences with respect to M85 OT2006-1 and M31 RV or only the result of a slower evolution 
is unconstrained. 
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Kennicutt et al. 20031 ). SINGS comprises IRAC observations of 75 galaxies within 30Mpc 



(mean distance 9.5 Mpc) with sensitivities similar to the observations presented here. We 
predict that a search in the archival data using the described color criterion might lead to 
the detection of as much as 1-10 new LRN events. A similar number of events could be 
found in a large survey of nearby galaxies using the deep imaging mode of Akari (ASTRO-F; 



Matsuhara et al. 20051 1 . 



We close the paper by emphasizing that the search for stellar mergers can be expected 
to have a substantial impact over the next years. Similar to the history of the field of 
gamma-ray bursts, a large increase of well studied events is anticipated to lead to exciting 
new insights into this emerging family of enigmatic transients. 
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Table 1: Log of observations. 



Start UTC 


Instrument 


Filter 


Exposure 


Flux density 


Brightness 






H 


[s] 




[Vega mag] 


2006 July 6.28 


Keck/NIRC 


2.1 


1200 


13±1 


19.3±0.1 


2006 July 7.06 


Spitzer/IRAC 


3.6 


3000 


35±6 


17.2±0.1 


2006 July 7.06 


Spitzer/IRAC 


4.5 


3000 


40±5 


16.6±0.1 


2006 July 7.06 


Spitzer/IRAC 


5.8 


3000 


46±6 


16.0±0.1 


2006 July 7.06 


Spitzer/IRAC 


8.0 


3000 


40±5 


15.5±0.1 


2006 July 2.12 


Spitzer/IRS 


15.6 


3000 


20±7 


14.9185 


2006 July 2.12 


Spitzer/IRS 


22.0 


3000 


25±20 


14.0iJ;« 



Table 2: Inferred black body parameters. 



Source 


-^peak 


T 

- 1 - eft , peak 




~^late a 


^eff.latc 


-Rlate 


Rl&te/t 




[xlO 5 L ] 


[xlO 3 K] 


[xlO 3 R ] 


[xlO 5 L ] 


[xlO 3 K] 


[xlO 3 R ] 


[km s _1 ] 


M85OT2006-l b 


~ 50 


~ 4.6 


~ 3.6 


2 9 +0 - 4 


0.95 ±0.15 


20±!j 




M31 RV C 


~ 8 


rvj 4 


~ 2 


~ 0.6 


~ 1 


~ 8 


~ 920 



a at t ~ 180 days for M85 OT2006-1 and t ~ 70 days for M31 RV 

b peak values from ( Kulkarni et al 20061 ). Late time values this paper 

°peak luminosity from ( Rich et al,1989( ). Remaining values from ( Mould et al. 199oh 



Fig. 1. — left: IRAC color-composite image of M85 using 3.6, 5.8 & 8/xm photometry 
obtained ~ 6 months after the eruption. The position of M85 OT2006-1 is indicated by the 
circle, right: Zoom into the location of the transient after subtraction of the galaxy light. 
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100 

JD-2453000 [days] 

Fig. 2. — if-band light curve, t = corresponds to the discovery date, 2006 January 7. 
Data at t < 120 days were reported in (IKulkarni et al 20061 ). The early-time light curve can 
be approximated by a linear decay with a slope of ~0.02mag per day. The late-time NIRC 
photometry suggests a significant flattening between 120 and 180 days. 



2 



4 6 8 io 

wavelength (/xm) 



20 



Fig. 3. — Near and mid- infrared flux densities observed in July 2006 (crosses). The solid 
line shows a black body with T e g = 1030 K fitted to to 2.1— 22/zm photometry. Omitting 
the K-band suggests a lower T eff (900 K; dashed line) and an additional hotter component. 
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[5.B]-[8.0] K s -[3.6] 



Fig. 4.— left: Color-color diagram of [3.6]-[4.5] vs. [5.8] — [8] for M85 OT2006-1 (square), 
point-like sources in the S-COSMOS fiel d fdots; ISanders et al. 20071 . courtesy of the S- 
COSMOS team) and T-dwarfs (triangl es; IPatten et al. 2006h . The locations of stars and 



AGN (contours from IStern et al. 20051 ) are indicated. Unresolved galaxies constitute the 
distribution of sources with [5.8] — [8]>1 mag and [3.6] — [4.5]<0.5 mag. Only T-dwarfs have 
similar colors as M85 OT2006-1 after ~ 180 days. Black body colors follow the solid line, 
right: Same for [3.6] — [4.5] vs. K s — [3.6]. Strong near-infrared absorption bands in T-dwarfs 
lead to a larger discrepancy compared to a pure black body spectral shape. 



